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Abstract 

A multiple input multiple output (MIMO) two-way relay channel is considered, where two sources 
want to exchange messages with each other using multiple relay nodes, and both the sources and relay 
nodes are equipped with multiple antennas. Both the sources are assumed to have equal number of 
antennas and have perfect channel state information (CSI) for all the channels of the MIMO two-way 
relay channel, whereas, each relay node is either assumed to have CSI for its transmit and receive channel 
(the coherent case) or no CSI for any of the channels (the non-coherent case). The main results in this 
paper are on the scaling behavior of the capacity region of the MIMO two-way relay channel with 
increasing number of relay nodes. In the coherent case, the capacity region of the MIMO two-way relay 
channel is shown to scale linearly with the number of antennas at source nodes and logarithmically with 
the number of relay nodes. In the non-coherent case, the capacity region is shown to scale linearly with 
the number of antennas at the source nodes and logarithmically with the signal to noise ratio. 

I. Introduction 

Relay channels are the most basic building block for cooperative and multihop communication in 
wireless networks. In a relay channel, one or more nodes, without data of their own to transmit, help a 
source destination pair communicate. The origins of the relay channel - as a three terminal communication 
channel - go back to Van der Meulen [1]. Despite the passage of time, the capacity of even the most 
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basic relay channels is still unknown. Nonetheless, bounds derived in [1], [2] show that using a relay, it 
is possible to increase the reliable rate of data transfer between the source and the destination. 

Motivated by the capacity improvements obtained by using multiple antennas at the source and the 
destination for point-to-point channels [29], recently, there has been a significant research focus on 
finding the capacity of the multiple input multiple output (MIMO) relay channel, where the source, the 
destination, and the relay may have multiple antennas [3], [10], [11]. The capacity of the MIMO relay 
channel was first studied in [3], [13], where upper and lower bounds on the capacity of the MIMO relay 
channel are derived for the deterministic and the Gaussian fading channel. Improved lower bounds for 
the MIMO relay channel with Gaussian fading channel were provided by [10], where message splitting 
and superposition coding are used at the transmitter to improve the bounds provided in [3]. In [3], [10] 
only full-duplex relays (can transmit and receive at the same time) were considered. Upper and lower 
bounds on the capacity for the more practical Gaussian MIMO relay channel with half-duplex relays, 
where the relays cannot transmit and receive at the same time, were developed in [11]. The bounds in 
[3], [10], [11] indicate that with relays there is a potential capacity gain to be leveraged by using multiple 
antennas. 

In [l]-[3], [10], [11] only a single source destination pair is considered with a single relay node. 
For a practical wireless network setting, where there are multiple source destination pairs, the concept 
of cooperative communication has been recently proposed [4]-[7], where different users in the network 
cooperate by taking turns relaying each others data. Thanks to the spatial separation between users, 
cooperation between users provides a means to obtain and exploit spatial diversity gain, called cooperative 
diversity gain, which increases the achievable data rate between each source and its destination. Several 
different protocols have been proposed to exploit the cooperative diversity gain, e.g. amplify and forward 
(AF) [4]-[7], [13], decode and forward (DF) [15], [18], with half-duplex [14], and full-duplex assumptions 
[16]. 

Prior work on the relay channel mostly considers one-way communication, i.e. a source wants to 
send data to a destination. In most networks, however, the destination also has some data to send to the 
source, e.g. packet acknowledgements from the destination to the source, downlink and uplink in cellular 
networks. Consequently, there has been interest in the two-way relay channel, where the bidirectional 
nature of communication is taken into account [21]-[24]. The two-way relay channel was studied in [22], 
where upper and lower bounds on the capacity region were derived for a general discrete memoryless 
channel. 

The MIMO two-way relay channel was introduced in [21], where two terminals T\ and T2 want to 
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Fig. 1. Two way relaying Protocol 

exchange information with each other through a single relay node as shown in Fig. [1] and both T\ , T2 
and the relay node is equipped with multiple antennas. It was assumed in [21] that each node can only 
work in half-duplex mode and there is no direct path between T\ and T2. The communication protocol 
proposed in [21] for the MIMO two-way relay channel is as follows. In the first time slot, both T\ and 
T2 transmit simultaneously and the relay node receives the superposition of the signals transmitted by T\ 
and T2. In the next time slot, the relay node transmits an amplified version of the signal, received in the 
last time slot, to both T\ and T%, subject to a power constraint. Since both T\ and T2 know what they 
transmitted in the last time slot, both can remove the effect of their own signal from the received signal, 
to decode the other terminal's message. Thus, the MIMO two-way relay channel facilitates simultaneous 
communication between T\ and T2 without creating any self interference. This idea is reminiscent of 
network coding [28], though note that here the coding is done in analog domain rather than in digital 
domain. The MIMO two-way relay channel is also known by several other names in the literature, namely, 
bidirectional MIMO relay channel [24] and is also a special case of analog network coding [28]. 

In prior work, achievable rate region (region enclosed by the rates achievable on the T\ — > T2 and 
T2 — ► T\ links, simultaneously) expressions were derived for the Gaussian half-duplex MIMO two-way 
relay channel (fading coefficients as well as additive noise is Gaussian distributed) using AF [21] and 
DF [23], [24] at the relay node. A main conclusion derived in prior work [21], [23], [24], is that it is 
possible to remove the | rate loss factor in spectral efficiency due to the half-duplex assumption on the 
relay node. To the best of our knowledge, none of the achievable rate region expressions for the MIMO 
two-way relay channel meet the best known upper bounds [25] and therefore the capacity region of the 
MIMO two-way relay channel is unknown. 

In this paper we consider a MIMO two-way relay channel with multiple relay nodes, deriving upper 
and lower bounds on its capacity region with different channel state information (CSI) assumptions. We 
show that the upper and lower bounds are only a constant term away, as the number of relays K grows 
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large, K — > oo with probability 1. Thus, we characterize the scaling behavior of the capacity region of 
the MIMO two-way relay channel as the number of relay nodes grow large. Our approach is similar to 
the asymptotic (in the number of relays) capacity formulation of [12], [19], [20]. 

Our system model and the key assumptions are as follows. We assume that two terminals T\ and T2 
want to communicate with each other via K relay nodes. None of the relays have any data of their own 
and only facilitate communication between T\ and Ti. Both T\ and T2 are equipped with M antennas, 
while all the K relays have N antennas each. We consider a two-phase communication protocol, where 
in the any given time slot, for the first a, a e [0, 1] fraction of the time slot, both T\ and T2 transmit 
simultaneously and all the relays receive. In the rest 1 — a fraction of the time slot, all the relays 
simultaneously transmit and both T\ and T2 receive the signal transmitted by all relays. We assume that 
there is no direct path between T\ and T2 and that Ti, T2 and all the nodes (Ti, T2 and all relay nodes) 
can only operate in half-duplex mode. No direct path assumption is reasonable for the case when relay 
nodes are used for coverage improvement and the signal strength on the direct path is very weak. The 
half-duplex assumption is made since full-duplex nodes are difficult to realize in practice. We assume that 
both T\ and T2 have perfect CSI for all the channels of the MIMO two-way relay channel in the receive 
mode. This could be enabled through a combination of channel reciprocity and feedback, however, we 
do not explore the practicalities of this assumption in this work. We consider two different assumptions 
about the availability of CSI at each relay node. First we consider the case when each relay is assumed 
to have perfect CSI for its own transmit and receive channel states, which is denoted the coherent MIMO 
two-way relay channel. Second we consider the case where the relays are assumed to have no CSI for 
any of their channel states, which is denoted the non-coherent MIMO two-way relay channel. 

Under similar assumptions, capacity scaling results have been found in [19] and [20] for the one-way 
relay channel (when T2 has no data for Ti). With a single antenna at both Ti, T2 and each relay node, it 
is shown in [19] that the capacity of the one-way relay channel scales logarithmically in the number of 
relay nodes, as the number of relay nodes grow large. The capacity scaling result of [19] was extended in 
[20] to the case where the source and the destination are equipped with M antennas and the all the relay 
nodes are equipped with N antennas and it was shown that there is a M fold increase in the capacity 
compared to the single antenna nodes [19]. 

The main results in this paper are on the capacity scaling laws for the MIMO two-way relay channel. 
For the coherent MIMO two-way relay channel, the capacity region is given by the convex hull of 

Rl2 < ylogK + 0(l) 

M 

R21 < ylogK + 0(l) 
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with probability 1 as K — > oo, where R12 and R21 is the rate of information transfer from T\ — > T2 
and T2 — ► Ti, and we use the notation u(x) = Oivix)) if I^r4| remains bounded, as x — ► 00. For this 
result, the upper bound on the capacity region is obtained for all a G [0, 1] and an achievable strategy 
with a = \ is proposed to achieve the upper bound within a 0(1) term. 

For the non-coherent MIMO two-way relay channel, for a fixed a = \, i.e. T\ and T2 transmit and 
receive for same amount of time, the capacity region is given by the convex hull of 

M 

R12 < —\ogP R + 0{l) 
M 

R21 < y log + 0(1) 

with probability 1 as K — > 00, and Pr is the sum of the power available at each relay. 

The strategy we use for deriving the capacity region of the MIMO two-way relay channel as K — > 00 
is to obtain an upper bound using the cut-set bound [30] and then derive an achievable rate region that 
approaches the upper bound. For the coherent case, we propose the following achievable strategy. Both 
Ti and T2 transmit M independent data streams from their M antennas. Each relay node using its CSI, 
does match filtering for the channels experienced by the M data streams from Ti — > T2 and T2 — > Ti, 
simultaneously, and all the M streams from Ti are decoded jointly at T2 and vice versa. We show 
that this strategy achieves the capacity region upper bound within a 0(1) term without any cooperation 
between Ti and T2. For the non-coherent case, we propose an achievable strategy where both Ti and T2 
transmit M independent data streams from their M antennas. Since none of the relays have any CSI in 
this case, we propose an AF achievable strategy where each relay transmits a scaled version of received 
signal subject to its power constraint, similar to [20]. With this strategy, as K — > 00, the channel between 
Ti — > T2 and T2 — > Ti converges to an M x M matrix with independent and identically distributed 
entries that are Gaussian distributed and we show that the achievable rate region provided by this AF 
strategy is within a 0(1) term of the upper bound in the high signal-to-noise ratio (SNR) regime. 

From an analytical perspective our work is closely related to [20], which only deals with MIMO one- 
way relay channel. We summarize the key differences and improvements of the proposed work compared 
to the MIMO one-way relay channel capacity scaling result of [20] as follows. 

• We assume a sum power constraint across all the relays, which is a generalization of the individual 
power constraint considered in [20]. An individual power constraint might seem more reasonable 
from a practical point of view. We show that even with a sum power constraint, however, the upper 
bound on the capacity region of both the coherent and non-coherent MIMO two-way relay channel 
can be achieved by allocating equal power to all relay nodes. Thus, with a sum power constraint, as 
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K —> oo, the total power transmitted by all relay nodes remains bounded as opposed to [20], where 
it is unbounded. The optimal power allocation is similar to [20] from a practical perspective, since 
each relay node is required to transmit the same amount of power. 

• We upper bound the capacity of the coherent MIMO two-way relay channel over all possible two- 
phase protocols, i.e. over arbitrary a, while in [20] an upper bound is derived only for a = \. 

• Our achievable AF strategy for the coherent MIMO two-way relay channel allows all the relays 
to help all the data streams going from T\ and T2 and T2 to T\ as opposed to [20] where only 
K/M relays are allowed to help each data stream. Moreover, in our AF strategy joint decoding is 
performed at both the receivers in contrast to [20], where each data stream is decoded by a single 
receive antenna treating all other streams as interference. Due to both these advantages, our AF 
strategy provides with better achievable rate regions compared to [20] for any finite K and a better 
0(1) term as K — > 00. 

• For the non-coherent MIMO two-way relay channel, we derive an upper and lower bound for a = \, 
which differs by only a constant term at high SNR, while in [20] only an achievable AF strategy is 
provided without any upper bound for a = \. 

Our results show that with the MIMO two-way relay channel there is a improvement in the capacity 
scaling by a factor of 2, compared to MIMO one-way relay channel [20], for both the coherent and the non- 
coherent case. We show that with the MIMO two-way relay channel, both T\ and T2 can simultaneously 
communicate with each other at a rate which is equal to the maximum rate at which T\ can communicate 
to T2 if T2 was silent. Therefore as K — > 00, the MIMO two-way relay channel is shown to create two 
interference free parallel channels, one for T\ — > T2 and another for T2 — * T\, where on each channel 
a rate given by the maximum possible rate at which T\ can communicate to T2 link if T2 was silent 
(one-way communication [20]) is achievable. 

Organization: The rest of the paper is organized as follows. In Section [nj we describe the MIMO two- 
way relay channel system model, the protocol under consideration and the key assumptions. In Section 
[ill} we derive an upper bound on the capacity of the coherent MIMO two-way relay channel. In Section 
Jvj by using a simple combining operation at the relays, we derive the asymptotic achievable rate region 
for the coherent MIMO two-way relay channel and show that it is possible to achieve the upper bound 
on the capacity region of the coherent MIMO two-way relay channel within a 0(1) term. Section [v] 
summarizes and discusses the implication of the coherent MIMO two-way relay channel capacity region. 



For the non-coherent MIMO two-way relay channel, in Section VI-A we derive an upper bound on the 



achievable rate region. Section VI-B gives a result on asymptotic achievable rate region for the non- 
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coherent MIMO two-way relay channel using AF strategy at relays. We draw some final conclusions in 
Section IVlTl 

Notation: The following notation is used in this paper. The superscripts T * represent the transpose 
and transpose conjugate. M denotes a matrix, m a vector and rrtj the i th element of m. For a matrix 
M = [mi rri2 . . . m n ] by vec(M) we mean [raj . . . rn^] T . det(M.) and ir(M) denotes the 
determinant and trace of matrix A, respectively. K x (f(x)) denotes the expectation of function / with 
respect to x. || • || denotes the usual Euclidean norm of a vector. \ m is a m x m identity matrix. \X\ 



u(x) 
v(x) 



remains bounded, 



is the cardinality of set X. We use the usual notation for u(x) = 0(v(x)) if 
as x — > oo. A circularly symmetric complex Gaussian random variable with zero mean and variance 
a is denoted by x ~ CJ\f(0,a) and x\y ~ CM(0,a) denotes that given y, x is a circularly symmetric 
complex Gaussian random variable with zero mean and variance a. The variance of a random variable a 
is denoted by var(a). C MN denotes the set of M x N matrices with complex entries. x n w ' p ' 1 > y denotes 
that the sequence of random variables x n converge to a random variable y with probability 1. We use 
o = b to denote equality with probability 1 i.e. Prob.(a = b) = 1 and < is defined similarly. I(x; y) 

w.p.l ' " w.p.l 

denotes the mutual information between x and y and h(x) the differential entropy of x [30]. To define 
a variable we use the symbol :=. 

II. System and Channel Model 

In this section we describe the MIMO two-way relay channel communication protocol, followed by 
signal and channel models. Consider a wireless network where there are two terminals T\ and T2 who 
want to exchange information via K relays, as shown in Fig. [2j The K relays do not have any data of 
their own and only help T\ and T2 communicate. We assume that there is no direct path between T\ 
and Ti and that they can communicate only through the K relays. This is a realistic assumption when 
relaying is used for coverage improvement in cellular systems, since at the cell edge the signal to noise 
ratio is extremely low for the direct path. In ad-hoc networks, this occurs when two terminals want to 
communicate, but are out of each other's transmission range. 

We assume that both the terminals T\ and T2 have M antennas while all the K relays each have 
antennas. The terminals T±,T2 and all the relays operate in half-duplex mode i.e. cannot transmit and 
receive at the same time. The communication protocol is summarized as follows [21]. In any given time 
slot, for the first a fraction of time, called the transmit phase, both T% and T2 are scheduled to transmit 
and all the relays receive a superposition of the signals transmitted from T\ and T2. In the rest (1 — a) 
fraction of the time slot, called the receive phase, all the relays are scheduled to transmit simultaneously 
and both the terminals receive. 
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Fig. 2. Two way communication in two hops 




Fig. 3. Channel Model 



A. Channel and Signal Model 

In this paper we assume that all the channels are frequency flat slow fading block fading channels, 
where in a block of time duration T c (called the coherence time), the channel coefficients remain constant 
and change independently from block to block. We assume that T c is more that the duration of the time 
slot used by T\ and T2 to communicate with each other as described before. As shown in Fig. [3] let the 
forward channel between T\ and the k th relay be = [hi^ h.2k ■ ■ • hjv/fc] and the backward channel 
between k th relay and T\ be = [hlv h[y • • • hjfcZ-]. Similarly let the forward channel between the 
k th relay and T2 be Gf. = [gki gk2 ■ ■ • gfcAf] and the backward channel between T2 and the k th relay be 
Gl r) = [g$ gfk ■ ■ ■ SM k ]- We assume that H ^ G l r) e C NxM ,H^,G k G C MxN with independent 
and identically distributed (i.i.d.) CN(0, 1) entries to keep the analysis simple and tractable. The ideas 
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presented in this paper, however, apply to a broad class of channel distributions. 
In the transmit phase, the N x 1 received signal at the k th relay is given by 

rk = V ~M~ HfcX + V ~m~ G * + nfe (1) 

where x and u are the Mxl signals transmitted from Ti and T2 to be decoded at T2 and T\ respectively, 
with E{x*x} = E{u*u} = M, P is the power transmitted by T\ and T 2 and E fe and F k are the path 
loss and shadowing parameters from T\ and T2 to the relay, respectively. The noise n^, is a spatio- 
temporal white complex Gaussian random vector independent across relays, with E(nfcn£) = ct 2 Itv- 
Relay k processes its incoming signal to transmit aJVxl signal ^/jkt k (with E{t£tfc} = 1) in the 
receive phase so that the transmitted power is 7^. We assume a power constraint of P at both Ti and 
T2 and a sum power constraint of Pr across all the relays, i.e. (J2k=i Ik < -Pr). The Mxl received 
signal v and y at terminal T\ and T2 respectively in the receive phase, are given by 

K 

v = ^ V / 7^04 r) t fc +w (2) 
k=i 
K 

y = ^2 VlkPkG k t k + z (3) 
k=l 

where 7^ is the power transmitted by the k th relay, Q k , P k are the path loss and shadowing parameters 
from the k th relay to T\ and T2, respectively, while w and z are M x 1 spatio-temporal white complex 
Gaussian noise vectors with E(ww*) = E(zz*) = g 2 \m- 

The path loss and shadowing effect parameters E k ,P k , F k and Q k V k for the link between T\ — > T 2 
and T2 — > Ti, are assumed to be independent and identically distributed (i.i.d.) random variables, strictly 
positive, bounded and remain constant over the entire time period of interest. 

(r) (r)t 

Throughout this paper we assume that both T\ and T2 perfectly know {H^, Hjj. , G k , G k '} V k, 
k = 1, 2, , . . . K in the receive mode. To be precise, in the receive phase (i.e. when T\ and T2 receive 
signal from all the relays), T x and T 2 both know {H fc , G k } and {H^ r) , Gj^} V k, k = 1, 2, , . . . K. We 
also assume that no transmit CSI is available at T\ and Ti, i.e. in the transmit phase T\ and T 2 have no 
information about what the realization of and G k is going to be when it transmits its signal to all 
the relays in the transmit phase, respectively. 

In this paper we consider two different assumptions about the CSI at the relays. The first case we 
consider is the coherent MIMO two-way relay channel, where all the relays have CSI in the transmit as 
well as the receive phases. For the coherent MIMO two-way relay channel, in the transmit phase the k th 

(r) (r) 

relay knows the realization of H^, Gjj. and in the receive phase it knows the realization of Gk, Hjj. , 
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K Relay Terminals 



Fig. 4. Broadcast Cut 

which could be achieved through channel reciprocity or feedback. We also consider the non-coherent 
MIMO two-way relay channel where we assume that no CSI is available at any relay. 

III. Upper Bound on The Capacity Region of The Coherent MIMO Two-Way Relay 

Channel 

The main result in this section is an upper bound on the rate R±2 and R21 of reliable transmission 
from T\ to T2 and from T2 to T\, given by the next Theorem. 

Theorem 1: The capacity region of the coherent MIMO two-way relay channel is upper bounded by 

M 

lim R 12 < — log if + 0(1) 

K^oo w.p.l 2 

M 

lim R21 < — log if + 0(1). 

K^oo w.p.l 2 

Outline of the Proof: We start by first separating T\ and then T2 from the network and apply the cut-set 
bound [30] to upper bound the rate of information transfer between T\ — > T2 and T2 — ► Ti, respectively. 
Using the cut-set bound, we first show that the maximum rate of information transfer from T\ — ► T2 
(T2 — > Ti) is upper bounded by the maximum rate of information transfer between T\ (T2) and relays 1 
to K (the broadcast cut) and also by the maximum rate of information transfer between relays 1 to K 
and T2 (Ti) (the multiple access cut), Fig. |4]and Fig. [5] Then we use the capacity result from Section 
4.1 [29] to upper bound the maximum rate through the broadcast cut for the case when CSI is only 
available at the receiver (all relays) and all the relays collaborate to decode the information. Similarly, 
for the multiple access cut as shown in Fig. |5J we upper bound the maximum rate at which all relays 
can communicate to T2 and Ti by using the capacity result from Section 3.1 [29], when CSI is known 
both at all the relays and at T% and T2 and all the relays collaborate to transmit the information. 
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Remark 1: For the broadcast cut, the upper bound on the capacity of the MIMO one-way relay channel 
can be found in [20] which also trivially serves as an upper bound on R\2 and i?2i- It is easy to identify, 
however, that there is a gap in the proof of Theorem 1 [20]. In Theorem 1 [20], it is argued that 

J(s;n,r 2 , • • . ,rjf|ti,t2, . . . ,t K ) = I (s; ri,r 2 , ...,r K ) 

where tfc is a function of r^, k = 1, 2, . . . , K , to which a counterexample can be easily found. Thus, we 
do not use result of [20] directly and attempt a different proof which is quite similar to the one given in 
[20], but closes the gap. 
The formal proof is as follows. 

(r) (r) 

Proof: Throughout this proof we assume that both T\ and T 2 perfectly know Hfc,HL ,Gfc, GL , 
Ek, Ffr, Qk, and Pk for k = 1, 2, . . . K, in the receive phase, and the k th relay knows H^, G^, E^, Fk 

(r) 

in the transmit phase and H^. ,Gk,Qk, Pk hi the receive phase. For notational simplicity, we do not 

(r) (r) 

include H^, G) ,Ek,Fk, H k , Gk, Qk, Pk in the mutual information expressions. We clearly point out, 
though, whenever their knowledge is used to derive the upper bound. 

Broadcast cut - To prove the upper bound we make use of the cut-set bound (Section 14.10 [30]). 
Separating the terminal T\ from the rest of the network and applying the cut-set bound on the broadcast 
cut as shown in Fig. |4j 

R12 < I(x;ri,r 2 , . . . ,r x ,y|ti,t 2 , . ..t K ,u). (4) 
Applying the cut-set bound while separating the terminal T 2 , 

#21 < I(u;ri,r 2 ,...,r^,v|ti,t 2 ,...t^,x) (5) 
for some joint distribution p(x, ti, t 2 , . . . , t^, u). By definition of mutual information [30] 
J(x;ri,r 2 , . . . , r K , y |t 1; t 2 , ...t K ,u) = 7(x;ri,r 2) . . . ,r^|ti,t 2 , . . .t K ,u) 

+ J(x;y|ri,r 2 , ...,r K , ti,t 2 , . . .t^,u). 
Expanding the mutual information in terms of differential entropy, 

I(x;ri,r 2 ,...,rx|ti,t 2 ,...t A ',u) = ft(x|t 1 ,t 2 , . • .t^u) 

- /i(x|ri,r 2 , . . . ,r^,ti,t 2 , ...t K ,u). 

Since conditioning can only reduce entropy [30], 

I(x;ri,r 2 ,...,rx|ti,t 2 ,...t^,u) < h(x\u) 

- /i(x|ri,r 2 , . . . ,r\K,ti,t 2 , ...t K ,u). 
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Note that ti,t 2 , • • • , tj<- is a function of ri, r%, . . . , yk, which implies 
/(x;ri,r 2 ,...,rx|ti,t 2 ,...t jft ',u) < /i(x|u) 

- fr(x|ri,r 2 , . . . ,r K ,u) 

and hence 

J(x;ri,r 2 , . . . ,r^|ti,t 2 , . . .t^u) < I(x;ri,r 2 , . . . ,rjf|u). 
From ([3]), with knowledge of Pk and G&, V k at terminal T 2 , 

J(x;y|ri,r 2 , . . . ,r K , ti,t 2 , ...t K ,u) = J(x,z) 
where z is the AWGN noise. Since x and z are independent, /(x, z) = 0, and therefore 
I(x; ri, r 2 , . • • , rx, y|t 1 ,t 2 , . . . t K , u) < J(x; n, r 2 , . . . , r^|u). 

Note that 

I(x;ri,r 2 , . . . ,r K , \u) = / x; —=,—=,..., —=|u 



K VK VK 
h ( -^L, ^L, . . . , -^|u ) - /i ( ^L, -^L, . . . , ^L|x,u ) . (6) 

/r Vk \[k ) \vk Vk Vk 



Next we evaluate §6ty using ([I]). From ([I]), 



r fc = A / — H fe x + W — G^u + n*. 



(V) 

Now with knowledge of i 7 ^ and G! at each relay 



,fr 1 r 2 r K \ I PE X m /P£ 2 n 2 l PEk n K 

Since conditioning can only decrease entropy, 



,/ri r 2 PE! m /P# 2 „ n 2 *JC 

HTf'Tf' •••'7f |u J- /l U™ 1 x + 7f'V™ H2X + 7f'---'Vw H - x + 7f 

With perfect knowledge of Ffc and H^, G^ at each relay 

, , ri r 2 \ / ni n 2 

h — =, — =, — = x, u = ft 



>K VK VK J WK VK VK 
and using ([6]) it follows that 



ni n 2 



J A^TT J * * * J 



K VK VK 
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Using the capacity result from Section 4.1 [29] when CSI is only known at the receiver, the R.H.S. can 
be upper bounded by 



\ K k=l I 



and the maximum is achieved when x is circularly symmetric complex Gaussian with E(xx*) = Im- 
Thus, it follows that 

/ 1 K p \ 

/(x;ri,r 2 ,...,r x |u) < logdet I A/ + ^ — — H£H fe . (7) 

V K k=l J 

Similarly, by interchanging the roles of x and u and replacing Ek with F k and with G&, 

J(u; n, r 2 , . . . , rxlx) < log det (l M + ^ ^ ^G[ r) *G[ r) ) . (8) 

\ K k=l / 



Using the strong law of large numbers 

K 



k=l 

and 

K 



Hm I y ^H£,H, ^ ^E{£ fc H£H fc } 



hm 1 V ^* G < r) *G< r) ^ ^E i^GpG^ j . 

fc=i 

Since E{H*H fc } = E{Gj. r) *G^ r) } = AHE M and let E{E k } = E{F k } = \i, using Q, Q, Q and the 
fact that the sources Ti and T 2 transmit only for a fraction of the time in each time slot, it follows that 

lim R W < aMlog ( 1 + ) (9) 

and 

lim i? 21 < aMlog ( 1 + ) . (10) 

Since M , AT", P, n and <r 2 are finite integers, as K —> oo 

lim fl 12 < aMlog(if) + 0(1) (11) 

If— »ac w.p.l 

and 

lim ft 2 i < aMlog(iT) + 0(1). (12) 

K^oo w.p.l 

Multiple access cut - Again by using the cut-set bound, we bound the maximum rate of information 
transfer i?i 2 (R21) from T\ — ► T 2 (T 2 — > T\) by the maximum rate of information transfer across the 
multiple access cut as shown in Fig. [5] Using the cut-set bound, i?i 2 and R21 are bounded by 

R12 < J(x,ti,t 2 , . . . ,tfir;y|u) (13) 
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R 21 < J(u,t 1 ,t 2 ,...,tjr;v|x). (14) 

By definition of mutual information 

7(x,t 1 ,t 2 , . • • ,t K ]y\u) = h(y\u) - /i(y|ti,t 2 , . . . ,t K ,u) 

+ /i(y|ti,t 2 , . . . ,tjc,u) - /i(y|ti,t 2 , . . . ,tx,x,u). 
Note that given ti, t 2 , . . . , tx, y is independent of x and u, thus 

/i(y|ti,t 2 , . . . ,t K ,x,u) = /i(y|ti,t 2 , . . .,t K ,u) = /i(y|ti,t 2 , . . . ,t K ). 

Therefore 

I(x,ti,t 2 , . . . ,t K ;y\u) = h(y\u) - h(y\t 1 ,t 2 , . . -,t K ). 
Since conditioning can only reduce entropy, 

J(x,ti,t 2 , . . . ,tfi:;y|u) < h(y) - h(y\t 1 ,t 2 , . . . ,t K ) 

= I(ti,t 2 ,...,tftr;y). 



Hence from ( pL3| ), 

R12 </(ti,t 2) ...,t A -;y). (15) 

Following similar steps 

flai < /(ti,t 2 ,...,t x ;v). (16) 

Clearly R\ 2 , R 2 \ are bounded by the maximum rate of information across the multiple access cut Fig. [5] 
Next, we compute J(ti, t 2; . . . , t^; y). Recall From ([3]), that the received signal y at T 2 is 

K 



y = VikPkG k tk + z 



fc=i 
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Note that 

y 



7(ti,t 2 , . . .,t K ;y) = I ( ti,t 2 , • • • ,t K ; 



Dividing y by \/K, the scaled signal is 

K 



This can also be written as 

y i 



K VK L 



Fid VP 2 G 2 ... JP K G 



'K 



[\/7iti^/72~t 2 • • • Vik^kV + 



K 



Note that $ is a M x matrix. Now assuming that all the relays know G& Vfc (allowing cooperation 
among all relays), with total power available across all relays bounded by Pr, from Section 3.1 [29], 

min{NK,M} , (K\v\\ 

/^ti,t 2 ,...,t x ;-^J < max|0,log^^^j| (17) 



where A;, Z = 1,2,..., min {NK, M} are the eigen values of matrix and v is chosen such that 

min{NK,M} , ^ n 

max jo, ^ - — | = P R . 

By definition, = ^ J^feLi PfcGfcG^. From the strong law of large numbers 

1 K 

lim V P fc G fe G* fc ^ E {PfeGfeG^} = E {P k } E {G fe G£} = »NI M 

k=l 

since E{G fe G*} = NI M and /z :=E{P k }. Therefore, it follows that 

\i = Nn Vt = l,2,...Af. 

which implies 

/Pr^ 1 

i/ = 

\M Nfj. 

and 

( y s 

/ ( ti,t 2 , . . . ,tjsr; 



^ (KNp(P R 1 



Since M, N, Pr, a 2 and \i are all finite, 

lim j(ti,t 2 ,...,tA-;-^=) < Mlogif + 0(1). 



Moreover, since the relays transmit only for (1 — a) fraction of time in any given time slot, 

lim R 12 < (1 - a)I ( ti,t 2 , ...,t K ; -?=) < (1 - a)MlogK + 0(1). (18) 

^-+00 yj.p.l \ VKJ 
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Similarly, we can derive a bound for R21 by using ([2]) and ( [T6| ), 



lim R21 < (1 - a)/ ti,t 2 , . . . ,t K ; -= < (1 - a)M\ogK + 0(1). (19) 



Combining (lit, (12), ( 18 ) and ( 19 > 



lim R12 < min{a, 1- a}M log K + 0(1) 

K^co w.p.l 

lim i?2i < min{a,l - a}M log K + O(l). 
K^oo w.p.l 



Since a G [0, 1], min{a, 1 — a} < h, therefore 



M 

lim i? 12 < — log if + 0(1) 

K ^00 ui.p.l 2 

M 

lim i? 2 i < —logiC + 0(1). 

K^oo w.p.l 2 



Discussion: In Theorem [TJ we obtained upper bounds on Ru and R21 by using cut-set bound on 
the broadcast cut (Fig. |4|) and the multiple access cut (Fig. [5]). For the broadcast cut, the upper bound 
corresponds to the case when the transmitter T\ or T2 has no CSI while all the relays collaborate to 
decode the message sent by T\ or T2 with perfect CSI, while the upper bound in the multiple access 
cut corresponds to the case when all the relays collaborate to transmit data to T\ or T2 using all their 
NK antennas with transmit CSI available at all relays. An important point to note is that the upper 
bound obtained in Theorem [T] is for any arbitrary a, which implies that the upper bound is valid for all 
two-phased MIMO two-way relay channel protocols and not for only a = 1/2 as is the case in [20]. 

In the next section we illustrate a simple amplify and forward (AF) strategy whose achievable rate is 
a constant term away from the upper bound. 

IV. Lower Bound on The Capacity Region of The Coherent MIMO Two- Way Relay 

Channel 

In this section we propose an AF strategy to achieve the upper bound obtained in Theorem [T] on the 
capacity region of the MIMO two-way relay channel within a constant term. The motivation to consider 
AF is because with DF, at each relay, to decode Ti's message T2's message is treated as interference 
and vice-versa, which implies that the achievable rate region with DF is same as that of the achievable 
rate region for the multiple access channel [30]. Since the achievable rate region of the multiple access 
channel is strictly less than the upper bound derived in Theorem [T] one cannot hope to achieve the upper 
bound given by Theorem [T] using DF protocol. With AF, however, each relay processes the received 
signal using its CSI and transmits it to T\ and T2 in the receive phase without any decoding. Since both 
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T\ and T 2 know what they transmit, (i.e. T\ knows x and T 2 knows u) with perfect receive CSI, both T\ 
and T 2 can cancel the contribution of their own transmitted signal from the received signal and decode 
other terminal's message without any self interference. 

Before discussing the MIMO two-way relay channel with multiple relays, let us first consider the case 
of a MIMO one-way relay channel (i.e. T 2 has no data for T\) with only one relay. For this case, the 
optimal AF strategy to maximize mutual information at the destination is to multiply V2DU1 to the 
signal at the relay, where the singular value decomposition of Hi is UiDiV^ and Gi is U2D2V2 and 
D is a diagonal matrix whose entries are chosen by waterfilling [31]. Finding the optimal AF strategy 
for the multiple relay case is a non-trivial problem and has not been found to the best of our knowledge. 
Moreover, the two-way nature of our problem makes it even more difficult to find the optimal AF strategy. 

To obtain a lower bound on the achievable rates for the MIMO two-way relay channel we propose 
a dual channel matching AF strategy in which relay k multiplies -4= ( G£,H£ + J3.^*G k ) to the 



received signal and forwards it to T\ and T2, where (3k is the normalization constant to satisfy the power 
constraint. In this AF strategy, each relay tries to match both the channels which the data streams from 
T\ to T2 and T2 to T\ experience. In dual channel matching, the complex conjugates of the channels 
are used directly rather than the unitary matrices from the SVD of the channels [31]. This modification 
makes it easier to compute the achievable rates for the MIMO two-way relay channel. 

Together with dual channel matching we restrict the signal transmitted from T\ and T2, x and u, 
respectively, to be circularly symmetric complex Gaussian with covariance matrix E{xx*} = E{uu*} = 
Q with tr (Q) = M (to meet the power constraint) to obtain an achievable rate region for the coherent 
MIMO two-way relay channel. Moreover, we use 01 = \ i.e. T\ and T2 transmit and receive for same 
amount of time. The achievable rates R\ 2 and R 2 \ using the above described AF strategy are given by 
the following Theorem. 

Theorem 2: For the coherent MIMO two-way relay channel, the achievable rates are given by 

lim R 12 = ^ log (if) + 0(1) 

K^oo w.p.l A 

lim R 21 = ^ log (if) + 0(1) 

A^oo w.p.l z 

with no cooperation required between T\ and T 2 . 

Proof: From M, the received signal at the k th relay is given by 



M k + V M 



rfc = \/^H fe x+ A /^G^u + n fe . (20) 
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Using dual channel matching, at relay k the transmitted signal t& is given by 



G*H* + H^'G^* 



where is to ensure that E{t£tfc} = 1. The received signal at T2 is given by 



K 



y = ^2 VikPkG k t k + z. 



Expanding (22i using (20 1 and (21 

K 



E 

fc=i 



hkP k PE k 

M h 



k=l 



A' 



hkP k PF k 



k=l 



Mf3 k 



(21) 



(22) 



(r)* r ~ 1 (r)* 



Gjb(GK + H^'G 



K 
fc=l 



'^G fc (GK + H?*Gf)n t+ z. 



B,, 



Since u and the channel coefficients H^, , G^, are known at T2, V A;, the second term can be 

removed from the received signal at T2. Moreover, as described before x is circularly symmetric complex 
Gaussian vector with covariance matrix Q, thus the achievable rate for T\ to T2 link is [29] 



R 



1 



12 



log det I M + 



AQA* 



Ejfc=l B fcB* + 1 M 



since Ejnfcn^,} = E{zz*} = Im, V k. Using Q = -Aj and dividing the numerator and denominator by 



R 



12 



log det I M + 



K 



A A* 

if A 



Note that as K — > 00, the contribution from 4>Im can be neglected and it follows that 



R12 = \ log det ( I M + 



K 



A A" 

if A 



Using equal power allocation 7fc = to satisfy the total power constraint of Pr across all relays 



I(x; y) = log det I M + 



K 



A A* 

K K 



where A = £f =1 



i\Pg fc 



G fc (G*H*+HpG^ 



M it Ef =1 B,B* 



Hi. and B* 



:G fc (G*H*+HpG^ 



Using the strong law of large numbers, 



and 



A w.p.l / P 2 T 

lim — ► \ T7 kM 1m 

K^oo K V M 



A* to.p.l / P 2 

lim — — > \ —km Im 

K^oo K V M 
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since E {G k G* k } = E {H£H fc } = M1 M V k, E {GfcGp} = E {HpH fc } = 0I M V k and P k , E k , [3 k 
are i.i.d. with k = E j J Pk ^ k j . Moreover, from the strong law of large numbers 

for some finite 6>, since B^B^. are i.i.d. for each A; and each entry of B& has finite variance. Thus, using 
these approximations, 

1 , / KPk 2 M 2 
hm R 12 = - log det I M H ^ Im 

K^oo w.p.l Z \ u 

,. „ Af. A KPk 2 M 2 

hm i?2i = — log 1 + 

_ftT— >oo u>.p.l 2 y f 

Since M, P, k and 9 are finite, as K — > oo, 

M 

lim i? 21 = — logiT + 0(l), 

il— >oo w.p.l A 

and similarly 

M 

lim R 21 = —logK + 0(l). 

K— >oo ui.p.l z 



Discussion: Theorem [2] shows that the achievable rate in each direction T% — > T 2 or T 2 — > T\ with 
the coherent MIMO two-way relay channel using dual channel matching at each relay is given by 
^ log (K) + 0(1) as K — > oo. More importantly, Theorem [2] also shows that both T\ and T 2 can 
simultaneously transmit at rate 4f log (K) + 0(1), without affecting each other's data rate and without 
requiring any cooperation between themselves. The result can be interpreted as follows. With dual channel 
matching, the transmitted signals from T% and T 2 are coherently added by all relays and the equivalent 
array gain of log K is obtained at both the receivers T\ and T 2 for all the M data streams transmitted 
from T 2 and T\. Moreover, with perfect channel knowledge, T\ and T 2 can cancel the self interference 
their own transmitted signals create, which enable T\ and T 2 to simultaneously achieve the rate of 
4^- log (K) + 0(1), without requiring any cooperation. 

Recall that for the MIMO one-way relay channel an AF strategy was proposed in [20] to achieve 
the upper bound within a constant term. With the AF strategy of [20], M independent data streams 
are transmitted from T\ and all the relays are divided into M sets with each set helping a particular 
data stream and independent decoding of data streams is employed at the receiver. Compared to the AF 
strategy of [20], with dual channel matching all relays participate in transmission of all data streams from 
T\ to T 2 and T 2 to T\ and thus provides a better achievable rate region. Moreover, joint decoding of data 
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Fig. 6. Asymptotic Capacity of the coherent MIMO two way relaying 



streams at respective receivers with dual channel matching removes the adverse effect of inter-stream 
interference which is caused due to independent decoding of different data streams in [20]. Thus it is 
clear that dual channel matching improves the achievable rate regions as compared to the AF strategy of 
[20]. 

V. Coherent MIMO Two-Way Relay Channel Capacity Region 



Combining the results from Section III and Section IV we establish the following characterization of 
the scaling behavior of the capacity region of the coherent MIMO two-way relay channel. 

Theorem 3: Neglecting the 0(1) term, the capacity region of the coherent MIMO two-way relay 
channel is given by the convex hull of 

lim R 12 = ^log(K) 
M 

lim R 2 i = -log(K) 

K^oo w.p.l A 

where R\ and R2 are the rate of information transfer between T\ — > T2 and T2 — > T\, respectively. 
Proof: Follows from Theorem [T] and [2] ■ 
Discussion: The capacity region of the coherent MIMO two-way relay channel is illustrated in Fig. 
[6] Combining Theorems [T] and [2j the sum capacity (sum of R12 and R21) of the coherent MIMO two- 
way relay channel is given by Mlog (K) + 0(1), as K — > 00, which is exactly double of the capacity 



21 



achievable in each direction T\ — > Ti or T2 — > Ti [20]. The 0(1) term in the upper and lower bound can 
in general be different and hence we characterize the exact capacity up to a 0(1) term. An important 
implication of Theorem [2] is that with enough relays, the dual channel matching strategy is optimal in 
the sense of achieving the right capacity scaling. Therefore, neglecting the 0(1) term, what this result 
shows is that with the coherent MIMO two-way relay channel, one can communicate at rate 4^ log K 
from T\ — > T2 while simultaneously communicating at rate M- log A" from T2 —*T\. 

Recall that we obtained the lower bound in Theorem [2] by fixing a = \ i.e. T\ and T2 transmit and 
receive for equal amount of time. Since this lower bound is only a 0(1) term away from the upper 
bound, allocating equal amount of time for the transmit and the receive phase is optimal for the coherent 
MIMO two-way relay channel. 

From Theorem [2j it is also clear that the upper bound on the capacity region of the MIMO two-way 
relay channel is achievable within a 0(1) term without any cooperation between T\ and T2. This is 
significant since the upper bound is for some joint encoding between T\ and T2. This is made possible 
because with channel knowledge, both T\ and T2 are able to cancel off the self interference. 

Compared to the asymptotic capacity result for MIMO one-way relay channel [20], our results show 
that with the coherent MIMO two-way relay channel one can remove the \ rate loss factor on the capacity, 
which comes from the half-duplex assumption on the terminals and relays. Therefore with the coherent 
MIMO two-way relay channel it is possible to can achieve unidirectional full-duplex performance with 
half-duplex terminals. 

To compute the capacity of the coherent MIMO two-way relay channel, we assumed that CSI was 
available at both T\ , T2 and each relay stage, which is a very strict requirement to meet in practice. To do 
dual channel matching for the coherent MIMO two-way relay channel, the k th relay needs to know the 

(r) (r) 

realization of and G) of the transmit phase and the realization of G^ and of the receive phase. 
To cancel the self-interference and to detect the incoming signal in the receive phase, both the terminals 
T\ and T2 need to know the realization of G&} V/c for the transmit phase and {H^ , Gl } Vfc for 
the receive phase. In practice, this is a very strict and challenging requirement, but by sending training 
sequence and using standard channel estimation techniques together with intelligent channel information 
feedback algorithms, all the nodes can learn the required receive channel coefficients with good enough 
accuracy. 

For example, by sending training sequences from T\ and T%, the k th relay can learn and G^ 

(r) 

in the transmit phase (when T\ and T2 transmit signals to all the relays). Learning the Gfe and 
realization at the k th relay for the receive phase (when all the relays transmit and both T\ and T2 receive) 
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is more challenging. In time-division duplex system, however, by employing calibration at transmitter and 
receiver, the forward and backward channel can be assumed to be reciprocal in which case the realization 
of ; and G k for the receive phase is approximately equal to the realization of H^. and G k for 
the transmit phase. Instead if a frequency-division duplex (FDD) system is used, assuming block fading 

(r) 

channel, Gk and H k can be learnt at each relay for receive phase, by feeding back the information 

(r) 

about Gfc and H k from T\ and T<± in the transmit phase, learnt in the last receive phase at T\ and T<i- 
To decode the incoming signal and to cancel the self interference, T\ and T2 needs to know the 

(r) (r) 

realization of Gk,Hk of the transmit phase and the realization of G k ,H k of the receive phase. 

By sending training sequences from all the relays to T\ and T2, T\ and T2 can learn the realization 

frl 

of H k ',G k of the receive phase, respectively. At the start of receive phase each relay knows the 

(r) (r) 

realization of Gfc,Hfc,G^, ,H k , therefore if each relay transmits quantized channel information about 

(r) (r) 

Gfc,Hfc, G k , Hi using strategies such as Grassmannian codebook [32] etc. to T\ and T2, both T\ and 
T2 can learn the required CSI in the receive phase. 

VI. Non-Coherent MIMO Two-Way Relay Channel 

In the last section we derived the scaling behavior of the capacity region of the MIMO two-way relay 
channel when both T\ and T2 have receive CSI while all the relays have perfect transmit and receive CSI. 
It is well known, however, that acquiring accurate CSI in a real-time communication is a challenging 
problem (large overhead and complexity) and guaranteeing near perfect CSI is almost impossible in 
practice. Therefore in this section we study the scaling behavior of the capacity region of the MIMO 
two-way relay channel when CSI is only available at T\ and T2 in the receive phase and no CSI is 
available at any of the relays. Furthermore, for this case we fix a = \, i.e. T\ and T2 transmit and 
receive for equal amount of time (transmit phase is equal to receive phase). 

For the non-coherent MIMO two-way relay channel, we first upper bound the achievable rates from 
T\ — » T2 and T2 — > T\ using the cut-set bound for the multiple access cut. Then using a simple AF 
strategy at each relay, we compute the achievable rates from T\ — > T2 and T2 — ► T\ which are shown 
to be within a 0(1) term from the upper bound in the high signal to noise (SNR) regime, thereby 
characterizing high SNR capacity. 

A. Upper Bound on The Capacity Region of The Non-Coherent Two-Way Relay Channel 

As proved in the section |ni| the rate of information transfer from T\ — > T2 (T2 — * Ti) is upper bounded 
by the rate of information transfer between all-relays put together and T2(T\) (multiple access cut). We 
evaluate this upper bound in the following Theorem, when CSI is not available at any of the relay. 



23 



Theorem 4: In the high-SNR regime (large Pr), the rates R\ 2 and R 2 \ from T\ — > T2 and T2 — > Ti 
for the non-coherent MIMO two-way relay channel are upper bounded by 

M 

lim R l2 < — log(P fl ) + 0(l) 

A^oo ui.p.l ^ 

M 

lim P 21 < — log(P R ) + C(l), 

K^co w.p.l £ 

where Pr is the total power constraint across all relays. 

Proof: Using the multiple access cut-set bound from ( p"5| ) and (16), we have 

R12 < I(ti,t 2 , • • • ,t K ;y) 

and 

R21 < I(t 1 ,t 2 , . . .,t K ;v), 

for some joint distribution p(ti,t2, . . . ,tx) and with no CSI at any relay. Recall from ^ that the 
received signal y is given by 

K 

k=i 

with power constraint Yl k =\ 7k < Pr- Using the capacity result from Section 4.1 [29] for no transmit 
CSI 

EQE* 

IP 

E = [v/pTGiv^Ga • • • v^G^] e C Mx7V/f 
and Q is the covariance matrix of 

[VTltiv^ta • • • V7M T 6 C*** 1 

with [ti t2 • • • t^] T circularly symmetric complex Gaussian and equivalent power constraint of tr(Q) < 
Pr and the maximum is achieved when Q = j^InkxNK- Therefore, using Q = -^InkxNK 

Pr 



I(ti,t 2 , ...,t K ;y)< logdet ( I M + 

where 



R 12 < logdet I M + NKg3 



k=l J 



From the strong law of large numbers, 

lim - V P k G k G% ^ E {P fc G fe G£} = E {P k } E {G fc G£} 



fc=i 



Since E{G fc G£} = AH[ M and let [i := E{P fc }, 



R12 < logdet fl M + — ^Im 

MJ.p.l \ C 
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Since \i and a 2 are finite, for large Pr 

Ru < M\ogP R + 0(l). 

w.p.l 

Since T\ and T2 transmit only for half the time (a = h) in any given time slot 

M 

R12 < — log P R + 0{1). 

w.p.l £ 



Similarly, it can be shown that 



M 

R21 < — logP R + 0(l). 

w.p.l A 



B. Lower Bound on The Capacity Region of The Non-coherent MIMO Two-Way Relay Channel 

In this subsection we compute achievable rates R12 and R21 for the non-coherent MIMO two-way 
relay channel using a simple AF strategy at each relay. The strategy is the following: with no CSI at 
any relay, each relay just normalizes the received signal to meet its power constraint and retransmits it 
in the receive phase. With CSI available at each destination T% (T2), self interference generated by T\ 
(T2) is removed from the received signal and the equivalent channel between T\ — > T2 (T2 — ► T\) for the 
non-coherent MIMO two-way relay channel is given by z~2k=i HfcGfc ^X^fc=i H^Gj^. As K — > 00, 
this channel is shown to behave as i.i.d. MIMO Gaussian channel. Then by using the capacity results 
from [29], we lower bound the capacity region of the non-coherent MIMO two-way relay channel. We 
show that with approximately same power used at T\(T2) and all relays (i.e. P Pr), the lower bound 
meets the upper bound in the high SNR regime (high P). 

The following Theorem gives the expressions for achievable R\2 and R21 pair, when each relay uses 
AF. 

Theorem 5: In the high SNR regime, the achievable rate region for the non-coherent MIMO two-way 
relay channel using AF strategy at each relay, is given by 

M 

lim R l2 = log(PR) + 0(l) 

K^oo w.p.l Z 

M 

lim R21 = — log(PR) + 0(l). 

K^oo w.p.l z 

Proof: Recall from ([T]) that the received signal at each relay is given by 



PEh / PFh „ ( r ) 



Vk = V lif HfcX + V k + nk (23) 

Therefore the average received signal plus noise power at each relay is given by N(P(Ek + Fk) + o -2 ). 
We assume that the k th relay knows the average received signal plus noise power N(P(Ek + Fj.) + a 2 ) 



25 



and transmits t k = f N(P(E k +F k )+<r 2 )) ) Vk to ensure tnat ^{^fctfc} = 1. With this normalization, from 
(|2]) and ([3]), the received signal at terminal T\ and T2 is given by v and y, respectively, where 



E\ IkQk Tj(r) . 

\ / — ; — -, ; sv-H; r k + w 
^ y N{P{E k + F k ) + a 2 ) k 



Substituting for from ([[} in the above equation 

K I 



y " NM(P(E k + F k ) + a^ ktik * 



K 



ST I lkPPkF k ( r) 

^ NM(P(E k + F k )+a^ k ^ k U 



+ g ^N(P(E k k + P F k ) + a^) Gknk + " 
Since T2 knows u and has perfect CSI, it can cancel the self interference. Removing the self interference 
from y and dividing both sides by yK, 

K 



, _ j_<T / IkPPkEk TT 

y " ^hx\ NM(P(E k + F k )+a^ ktikX 



/ IkPk n , 1 



Similarly Ti knows x and also has perfect CSI, therefore it can also remove the self interference. 
Removing the self interference from v and dividing both sides by yK 



„/ _ 1 V" / lkPQkF k H (r) r (r) ii 

NM(P(E k + F k )+a 2 ) k k 



1 K 
IF ^ 



IkQk „(r) , 1 



# V + i^) + k Vk 



As if — > 00, it can be shown that (Theorem 3 [20]) 
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\ k=i ^ 



IkPQkFk 



M(P(E k + F fc ) + a 2 ) 
IkPk 




{P{E k + F k ) + a 2 ) 



and 



to.p.l 1 | j k PP k E k \ 

> Kt^\M{P{E k + F k )+^)) 1 

w.p.l Iv^lB-/ IkPQkFk 1, 



where Aj^Cjj denotes row and column entry of A and C respectively and bj,dj denotes 
the element of b and d respectively, Ra = E{aa*} where a = vec(A) and Rc = E{cc*} where 
c = vec(C). 

This shows that the channel matrices A, C and the noise vectors b, d are i.i.d. Gaussian, therefore 
using results from Section 4.1 [29] with only receive CSI and no transmit CSI, the achievable rate R12 
(R21) of the Ti -► T 2 (T 2 -► I\) link for a = §, is given by 

Jim R 12 = l-E Uw {logdet (l M + S H » H »)} 

A— >oo ui.p.l Z I \ M /J 



Jim i? 21 = W llogdet (l M + ^H^)) 



A^oo w.p.'. 

where H w is an M x M matrix with i.i.d. CJ\f(0, 1) entries and 

X 2^fc=l ^ I (p( jBfc+Ffc ) + CT 2 ) j> 



Pi 



P2 



a- (ELi E { (p(eJ+S)+^) } + !) 
i_ ^a Tg r -y k PQ k F k 1 



A (EfcLl E { (P(E?+F*)+<r>) } + 
Note that pi and p 2 are effective SNRs. Denoting fj, = E{E k } = E{F k } = E{P k } = E{Q k } Vfe, and 

^ = E {p(P fc +A)+a 2 i = E {p(P fc +F fc )+^ Vfe ' 

P/j ^A 
^ Z^fe=i 7fe 



Pi = P2 



° 2 {}Xk=iik + i) 



Since the relay power is constrained by Ylk=i ^ k = 

= PP Rf i 
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Choosing P « Pr, 



Pi = P2 




since & 



•k,Fk,Pk,Qk- Vfc are bounded. Therefore 



K^oo K^oo w.p.l 



lim i?i2 = lim i?2i = Eh to \ logdet ( Im + 



Ma 2 



Pr 




In high SNR regime P m P R > 1, from [29], it follows that 



lim R\2 



lim i?2i 



■uj.p.l 



w.p.l 



-log(PjO + 0(l) 
|log(P R ) + 0(l). 



Discussion: In this section, we first obtained an upper bound on the capacity region of the non-coherent 
MIMO two-way relay channel using multiple access cut-set bound when CSI is only known at T\ and 
Ti. Then with the help of a simple AF strategy we provided a lower bound which is a 0(1) term 
away from the upper bound in the high SNR regime. We find that, contrary to the coherent case, with 
the non-coherent MIMO two-way relay channel, as the number of relay nodes grow large, the capacity 
region expression is independent of the number of relays and no coherent combining gain (array gain) 
is available when there is no CSI at any relay. Similar to the coherent case, however, it turns out that 
even in the non-coherent case both T\ and T2 can simultaneously transmit at a rate which is equal to 
the maximum possible rate at which they could have transmitted when there is no data flowing in the 
opposite direction. Therefore, the non-coherent MIMO two-way relay channel creates two orthogonal 
channels, one from T\ — ► T2 and another from T\ — ► T2 with rate 4f log Pr achievable on each link 
simultaneously, thereby removing the g rate loss factor because of half-duplex nodes. 

The lower bound provided by Theorem [5] shows that the achievable rate for the non-coherent MIMO 
two-way relay channel is same as the capacity of a point to point M x M i.i.d. Gaussian channel with 
receive SNR Pr, with perfect CSI at receiver and no CSI at transmitter and where i factor is due to 
the half-duplex requirement. This result is quite intuitive, since with absence of CSI at the relays, as 
K — > 00 the equivalent channel between T\ — ► T2 (T2 — ► 2i) converges to an M x M i.i.d. Gaussian 
channel and therefore the result follows from [29]. 

Compared to (Theorem 3 [20]), this result shows that with the non-coherent MIMO two-way relay 
channel it is possible to remove the | rate loss factor due to the half-duplex constraint and can achieve 
the same rate as promised by Theorem 3 [20] (for unidirectional communication), in each direction 
T\ — * T2 and T^—^Tx. This is again due to the fact that, with perfect CSI both T\ and T2 can cancel the 
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self interference terms their own transmitted signals generate and hence the received signal at T2 (Ti) 
when T2 is also sending information is equivalent to the received signal at T2 in [20], where there is 
no communication happening on T2 — > T\ link. Therefore there is a two-fold increase in achievable rate 
with the non-coherent MIMO two-way relay channel in comparison to [20]. 

VII. Conclusion 

In this paper we developed capacity scaling laws for the MIMO two-way relay channel under coherent 
and non-coherent assumptions. First we upper bounded the capacity region of the coherent MIMO two- 
way relay channel using the broadcast and multiple access cut-set bound. Then we proposed a dual 
channel matching strategy to obtain an achievable rate region for the coherent MIMO two-way relay 
channel. The achievable rate region was shown to be a 0(1) term away from the upper bound, as 
K —> 00. Hence we characterized the coherent MIMO two-way relay channel capacity region within a 
O(l) term as K — > 00. 

The dual channel matching strategy we proposed for the coherent MIMO two-way relay channel is a 
decentralized strategy, where each relay node does not cooperate with any other relay node and only uses 
its CSI to coherently match the channels which the streams from T\ and T2 experience. An interesting 
outcome of our analysis is that the dual channel matching strategy, which requires no cooperation between 
relays, achieves the capacity region upper bound which allows for full cooperation between relays, within 
a O(l) term. Thus, dual channel matching not only simplifies the practical protocol design, but also 
achieves capacity region upper bound within a 0(1) term. 

For the coherent MIMO two-way relay channel, there is a strict requirement that all the nodes need 
to know perfect CSI, which in practice can be quite challenging and resource consuming. Therefore we 
also considered the case when only T\ and T2 have perfect receive CSI and none of the relays have 
any CSI, which is referred to as the non-coherent MIMO two-way relay channel. For this case we upper 
bounded the capacity region using only the multiple access cut-set bound and fixing a= \ (i.e. Ti, T2 
and all the relays transmit for equal amount of time in each time slot). Then with the help of a simple 
AF strategy, we showed that in the high SNR regime the upper bound is achievable within a 0(1) term, 
and hence characterize high SNR capacity region of the non-coherent MIMO two-way relay channel. 
Thus, we showed that a very simple AF strategy which transmits the power normalized version of the 
received signal is an optimal strategy. The intuition behind this result is that by using AF with no CSI 
at any relay, the effective channel from the source to destination T\ — ► T2 or T2 — > T\ converges to an 
M x M i.i.d. MIMO Gaussian channel as K —>■ 00, which is similar to the effective channel considered 
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for the capacity region upper bound (using the multiple access cut). The upper and lower bound differ 
by 0(1) term because of the forwarded noise from each relay node. 

Compared to [19], [20], our capacity scaling results for the coherent and non-coherent MIMO two- 
way relay channel shows that with the MIMO two-way relay channel there is a two-fold increase in 
the capacity than unidirectional communication with large number of relays. Hence, the MIMO two-way 
relay channel helps in improving the spectral efficiency and unidirectional full-duplex performance while 
using half-duplex terminals. 
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